REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  Wo.  07010188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704  0188).  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION 

1  REPORT  DATE  IDD-MM-YYYY)  2  REPORT  TYPE 

13-01-2010  Journal  Article 

3.  DATES  COVERED  ( From  -  To) 

4.  TITLE  AND  SUBTITLE 

Application  of  a  Vanishing  Quasi-Sigma  Vertical  Coordinate  for  Simulation 
of  High  Speed  Deep  Currents  Over  the  Sigsbec  Escarpment  in  the  Gulf  of 
Mexico 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c  PROGRAM  ELEMENT  NUMBER 

060243 5N 

6,  AUTHOR(S) 

D.  Dukhovskoy,  S.  Morey,  Paul  Martm,  J.  O'Brien,  C.  Cooper 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNlt  NUMBER 

73-6057-08-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/J  A/7320-08-83 1 4 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy  St. 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

ONR 

Arlington,  VA  22217-5660 

11,  SPONSOR/MONITOR'S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


20100121308 


14.  ABSTRACT 

Kcccnt  observations  over  the  Sigsbeo  Escarpment  in  the  Gulf  of  Mexico  have  revealed  extremely  energetic  deep  currents  (near  1  m  s'1),  which  tire  trapped  along  the  escarpment.  Both 
scicmiric  imerest  and  engineering  needs  demand  dynamical  understanding  of  these  extreme  events,  and  can  benefit  from  a  numerical  model  designed  to  complement  observational  and 
theoretical  investigations  in  this  region  of  complicated  topography.  The  primary  objective  of  this  study  is  to  develop  a  modeling  methodology  capable  of  simulating  these  physical  processes 
and  apply  the  model  to  the  Sigshee  Escarpment  region.  Ihe  very  steep  slope  of  the  Sigsbec  Escarpment  (0.05  -0.1)  limits  the  application  of  ocean  models  with  traditional  terra  in- following 
(sigma)  vertical  coordinates,  which  may  represent  the  very  complicated  topography  in  the  region  adequately,  can  result  in  large  truncation  errors  during  calculation  of  the  hori/.onlal  pressure 
gradient.  A  new  vertical  coordinate  system,  termed  a  vanishing  quasi-sigma  coordinate,  is  implemented  in  the  Navy  Coastal  Ocean  Model  for  application  to  the  Sigshee  Escarpment  region. 
Vertical  coordinate  surfaces  for  this  grid  have  noticeably  gentler  slopes  than  a  traditional  sigma  grid,  while  still  following  the  terrain  near  the  ocean  bottom.  The  new  vertical  grid  is  tested 
with  a  suite  of  numerical  experiments  and  compared  to  a  classical  sigma-layer  model.  The  numerical  error  is  substantially  reduced  in  the  model  with  the  new  vertical  grid.  A  one-year, 
realistic,  numerical  simulation  is  performed  to  simulate  strong,  deep  currents  over  the  Escarpment  using  a  very-high-resolulion  nested  modeling  approach.  The  model  results  are  analyzed  to 
demonstrate  that  the  deep-ocean  currents  in  the  simulation  replicate  the  prominent  dynamical  features  of  the  observed  intense  currents  in  the  region. 


15.  SUBJECT  TERMS 

numerical  models,  deep  currents,  topographic  waves,  numerical  truncation  error,  vertical  discretization,  Sigsbec  Escarpment 


16.  SECURITY  CLASSIFICATION  OF' 


a  REPORT 

Unclassified 


b  ABSTRACT 

Unclassi  fied 


c.  THIS  PAGE 

Unclassified 


17.  LIMITATION  OF 
ABSTRACT 

UL 


18.  NUMBER 
OF 

PAGES 

16 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Paul  J.  Martin 


19b.  TELEPHONE  NUMBER  (Include  area  code ) 
228-688-5447 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std  Z39.18 


PUBLICATION  OR  PRESENTATION  RELEASE  REQUEST 


Pubkey:  5915  nrlinst5600.2 


•  REFERENCES  AND  ENCLOSURES 

Ref:  (a)  NRL  Instruction  5600.2 
(b)  NRL  Instruction  5510  40D 

Enel:  (1)  Two  copies  of  subject  paper 
(or  abstract) 

4.  AUTHOR 


(  )  Abstract  only,  published  ( 

(  ) Book  ( 

(  )  Conference  Proceedings  ( 

(refereed) 

(  )  Invited  speaker  ( 

(  X  )  Journal  article  (refereed)  ( 

(  )  Oral  Presentation,  published  ( 

(  )  Other,  explain 


)  Abstract  only,  not  published 
)  Book  chapter 
)  Conference  Proceedings 
(not  refereed) 

)  Multimedia  report 
)  Journal  article  (not  refereed) 

)  Oral  Presentation,  not  published 


3.  ADMINISTRATIVE  INFORMATION 


2.  TYPE  OF  PUBLICATION  OR  PRESENTATION 


STRN  NRUJ  A/7320-  OS-  8314 

Route  Sheet  No  7320/ 

Job  Order  No.  73-  6057-  08-  5 
Classification  X  U  C 

Sponsor  ONR 

approval  obtained  yes  X  no 


Title  of  Paper  or  Presentation 

Application  of  a  Vanishing  Quasi-  Sigma  Vertical  Coordinate  for  Simulation  of  High  Speed  Deep  Currents  Over  the  Sigsbee  Escarpment  in  the  Gulf  of 


Author(s)  Name(s)  (First,MI,Last),  Code,  Affiliation  if  not  NRL 

Dmitry  Dukhovskoy,  Steven  Morey,  Paul  J.  Martin,  James  J.  O’Brien,  Cortis  Cooper 


It  is  intended  to  offer  this  paper  to  the 


(Name  of  Conference) 


and/or  for  publication  in 


(Date,  Place  and  Classification  of  Conference) 
t/  L  <-  .  <  /  «£  /  /  /  f  / 

JourrraLof-Geapfoyeteat-RraeaFGh,  Uhclassified 

(Name  and  Classification  of  Publication)  (Name  of  Publisher) 

After  presentation  or  publication,  pertinent  publication/presentation  data  will  be  entered  in  the  publications  data  base,  in  accordance 
with  reference  (a). 

It  is  the  opinion  of  the  author  that  the  subject  paper  (is _ )  (is  not _ *)  classified,  in  accordance  with  reference  (b). 

This  paper  does  not  violate  any  disclosure  of  trade  secrets  or  suggestions  of  outside  individuals  or  concerns  which  have  been 

communicated  to  the  Laboratory  in  confidence.  This  paper  (does _ )  (does  not _ X)  contain  any  rfillitarily  critical  technology. 

This  subject  paper  (has _ )  (has  never  X  )  been  incorporated  in  an  official  NRL  Report  j  >  * 


Paul  J.  Martin,  7322 

Name  and  Code  (Principal  Author) 


"A 


(Signature) 


5.  ROUTING/APPROVAL 

C 

CODE 

, SIGNATURE 

DATE 

COMMENTS 

Author(s)  /,  A  .1  ,  - 

r  l  u  r  1  * ' 1 

~  7 

Ss  L.  -  <K 

C?/j7j  Ad  "  S' 

Need  bv  v  /  /  c  "  ^  ^  ^ 

Publicly  accessible  sources  used  for  this  publication 

Section  Head  / 

/  !  H  i  ,  <■*■ 

iiV1  _ 

Branch  Head 

Gregg  A.  Jacobs,  7320 

V//3 

Division  Head 

Ruth  H.  Preller,  7300 

Security,  Code 

1226 

I  1  O  6 

1 .  Release  of  this  paper  is  approved. 

2.  To  the  best  knowledge  of  this  Division,  the 

subject  matter  of  this  paper  (has _ _  ) 

(has  never  x  )  been  classified. 

1.  Paper  or  abstract  was  released.  ^ 

2.  A  copy  is  filed  in  this  office.  3  f7  ^  j 

Office  of  Counsel, Code 

1008.3 

i 

f  .■  ''A 

f 

ADOR/Director  NCST 

E,  R.  Franchi,  7000 

Public  Affairs  (Unclassified/ 

Unlimited  Only),  Code  703Q  4 

Division,  Code 

Author.  Code 

HQ-NRL  5511/6  (Rev.  12-98)  (e) 


THIS  FORM  CANCELS  AND  SUPERSEDES  ALL  PREVIOUS  VERSIONS 


Martin,  7322 


''(SignSfUre) 


.SIGNATURE 


Aj-vo^sr 


wp^sfsea 


VBpoVcha^ifir 
$C6ftf<frdh&  Proceedings 
;:.s  (fidtrtfeteed) 

)  Multimedia /sport 
)  Journal  "article  (hottsfenee.d) 

)  OrabPrestehtabon,  not  published 


•  ’■  *V,  ,  ■  (  r J..  ‘  ).,■  .  v' 

pal pfcr. 

'  y ^'(tf abstract) ‘ 


f$EhS)  :-  fl\  _■ 

V  thviW  speatfer, 

}  (rafere^tf) 

):Or*KRm&fentet!brw  published 


seen 


mm* g 


59 


INST  M30:  i 


jttoOrdGrNo.  73.CQ$7-0a-S  _ 

Classification  U 

Sponsor  ON*  c*?fyx(  (J  Q  * 
approval. obtained  yes  X  no 


..  *  .,.T,  ..  .  ■  - 

A ^ riaitd  for  SI mylaitton  of  High  Speed  Deep  Currents' Over  the  Stgsboo  Escarpment  In 


the  Gulf  of 


t> rft I R^av^tWSjf V en  Mdrey;  :PJra>J^fifljrttln,  Jame*x6*tif1eni  Cortis  Cooper 

;^:v::#0W^V'  ’"  '■  .;"'  '  .  r^.s 

■^>.: '■^■r.v. . .^:  ny^h^ys*ouritvP. 

‘  7  ')^f«pnro™a<fe  >n  ft/  '  ewe.'M/ 


M  fiXrco 


It  is  mtentieclio  offer  this  paper  to  the  _ _ 


(Nema  of  Conference) 


yj w 


-■?6 


(Dtte.PIbce  and  Glassification  of  Conference) 


r  ,..v,.  ,  .,  .  .  .^ .  (Narriebf  Publisher) 

Aft&r  ^^Utjfr^tii^Tipfef^^feTit  ^utilifcarfhiin/pres^rntatlon ^clata  will  be  entered  in  the  publications  data  base,  In  accordance 

with rafetiheeja)^’  .  ’  J,i,f:  7;;''. 

It  is  thfrojsfnfiM'SRIfe _ _ J(is  not  xl  classified,  In  accordance-iwith  reference  <b). 

This  papdrd^  hOT^ibiafe  artir  discIostjilp:6f  trade'Secretsi()r  suggestions  oiootside  individuals  or  concerns  which  Have  been 

comrTnjnicated'to^eiJtt»i«ory  Irr^oftfWence.  This paf>cr  (does  I fdoes  not _ X)  contain  any-rfllKartly  critical  technold&y 

This  subject  papdethas  _ _ )  (has  never  X  Tbedh incorporated  in  an  official  NRL  Report  y  *  " 


#7  Q*-i 


2*><r 


Pubfldy  accessible  source*  used  lor  Ihb  public®  ton 


- —  Personally  identifying  email  addresses  may  not 

_  _  be  posted  to  publicly  accessible  DON  web  sites 

per  SECNAV1NST  5720.47B,  end  (1>;  3. <3  (7) 


s^tnyrmstr 

'Brmrcri 


AT 


4-R-ca? 


■  ■ 

RirthffcPrellw!  TWp.v  '  ;.|v 

•  W  •  tjiff . .> 

■SScunty^CKe 

1226 

"DIm  o'  coureTCCSde' 

'"■  Wtffej,. 


?/ns  /t?5 


’iltirlO'S 


t.  Pa  par  or  absimtf  was  released 
2  „  A  copy  to  riled  H  this  office  J ' 


AOdR/bhector  NCST 

6.  R.  pranchli  TCOO 

.  _ ^  >■  :■  .1. 

Pubnc.AffelrsYUrtCfawrilfe*' 
Lfn/irrrited'OrtlyJi  Cod*  7jj3&.4 


DMsion,  Cod£- !  ■ 

_ ^  ~  ^ 

Aulhrai'CWS'/T'-  ;?:V;- 


■feaMum 


wfci  >  i  ^ito 


TrftSeaaa  apiHoiis 

2,  Toiho  b«*t  kno^iadgo  of  Tbte  Dtvri»oo,  me 

iiubjoct  manor  of  th«  papor  (^st _ 1 

(has  never ,  x  )  been  cMssrfitd. 


HQ-NRL  55 life titBV  ’f  JjsStfblr'- 

SS-Snd  C6.d3S 


THISTORM  CANCELS  i^lD  SuPHRSEDESXLl  PR6VIOWS  versions 


Ocean  Modelling  28  (2009)  250  -26S 


•A' 

fwm' 

I  I  SI-YII  K 


Contents  lists  available  at  ScienceDirect 

Ocean  Modelling 

journal  homepage:  www.elsevier.com/locate/ocemod 


Application  of  a  vanishing,  quasi-sigma,  vertical  coordinate  for  simulation  of 
high-speed,  deep  currents  over  the  Sigsbee  Escarpment  in  the  Gulf  of  Mexico 

Dmitry  S.  Dukhovskoy3  *,  Steven  L  Morey a,  Paul  J.  Martin  b.  James  J.  O'Brien a,  Cortis  Cooper c 

*  Center  far  Ocean  -  Atmospheric  Prediction  Studies,  The  Florida  State  University.  203S  E.  Paul  Dirac  Drive,  PM.  Jahnsan  8 Idg.  Suite  200,  Tallahassee.  FL  32306-2840,  USA 
b  Naval  Research  Laboratory,  Stennis  Space  Center,  MS  39S29,  USA 

c  Energy  Technology  Co..  Chevron .  6001  Bollinger  Conyon  Rd.  Rm  L-4240,  San  Romon,  CA  94583,  USA 


ARTICLE  INFO 


ABSTRACT 


Am'c/e  history: 

Received  2S  September  2008 
Received  in  revised  farm  26  January  2009 
Accepted  23  February  2009 
Available  anline  3  March  2009 


Keywords. 

Numerical  models 
Deep  currents 
Topographic  waves 
Numerical  truncation  error 
Vertical  discretization 
Sigsbee  Escarpment 
Gulf  of  Mexico 


Recent  observations  over  the  Sigsbee  Escarpment  in  the  Gulf  of  Mexico  have  revealed  extremely  ener¬ 
getic  deep  currents  (near  1ms  ,).  which  are  trapped  along  the  escarpment.  Both  scientific  interest 
and  engineering  needs  demand  dynamical  understanding  of  these  extreme  events,  and  can  benefit  from 
a  numerical  model  designed  to  complement  observational  and  theoretical  investigations  in  this  region  of 
complicated  topography.  The  primary  objective  of  this  study  is  to  develop  a  modeling  methodology  capa¬ 
ble  of  simulating  these  physical  processes  and  apply  the  model  to  the  Sigsbee  Escarpment  region.  The 
very  steep  slope  of  the  Sigsbee  Escarpment  (0.05-0.1)  limits  the  application  of  ocean  models  with  tradi¬ 
tional  terrain-following  (sigma)  vertical  coordinates,  which  may  represent  the  very  complicated  topogra¬ 
phy  in  the  region  adequately,  can  result  in  large  truncation  errors  during  calculation  of  the  horizontal 
pressure  gradient.  A  new  vertical  coordinate  system,  termed  a  vanishing  quasi-sigma  coordinate,  is 
implemented  in  the  Navy  Coastal  Ocean  Model  for  application  to  the  Sigsbee  Escarpment  region.  Vertical 
coordinate  surfaces  for  this  grid  have  noticeably  gentler  slopes  than  a  traditional  sigma  grid,  while  still 
following  the  terrain  near  the  ocean  bottom  The  new  vertical  grid  is  tested  with  a  suite  of  numerical 
experiments  and  compared  to  a  classical  sigma-layer  model.  The  numerical  error  is  substantially  reduced 
in  the  model  with  the  new  vertical  grid.  A  one-year,  realistic,  numerical  simulation  is  performed  to  sim¬ 
ulate  strong,  deep  currents  over  the  Escarpment  using  a  very-high-resolution  nested  modeling  approach. 
The  model  results  are  analyzed  to  demonstrate  that  the  deep-ocean  currents  in  the  simulation  replicate 
the  prominent  dynamical  features  of  the  observed  intense  currents  in  the  region. 

2009  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Observations  of  currents  in  the  deep  Gulf  of  Mexico  have  shown 
that  the  abyss  is  far  from  quiescent.  Strong  near-bottom  currents 
exceeding  1ms  1  have  been  measured  in  water  depths  of 
2000  m  over  the  Sigsbee  Escarpment  (SE),  south  of  New  Orleans 
(Hamilton  and  Lugo-Fernandez,  2001;  Hamilton  et  al.,  2003).  Evi¬ 
dence  of  extreme  deep  currents  in  the  region,  potentially  impor¬ 
tant  for  offshore  petroleum  exploration  and  extraction  activities, 
has  led  to  intensive  observational  efforts  in  the  area  (discussed  in 
Section  2)  and  has  sparked  interest  within  the  oceanographic  com¬ 
munity  (Hamilton.  1990;  Hamilton  and  Lugo-Fernandez,  2001 ;  Oey 
and  Lee,  2002;  Hamilton.  2007). 

Despite  a  growing  database  of  deepwater  current  measure¬ 
ments  and  analyses  of  the  observations,  mechanisms  generating 
the  deep  strong  currents  remain  unclear.  It  has  been  speculated 


•  Corresponding  author.  Tel.:  +1  8S0  644  1 168;  fax-  +1  850  644  4841 
E-mai7  addresses.  ddmitry#coaps.fsu,edu  (D.S.  Dukhavskay).  smorey^caaps. 
fsu  edu  (S.L  Marey),  paul.martin^nrlssc.navy.mil  (P.J.  Martin),  jabnen<3>caaps.fsu. 
edu  (j.J.  O'Brien),  cartcoopert^hevron.cam  (C  Cooper). 

1463-S003/S  -  see  frant  matter  ©  2009  Elsevier  Ltd.  All  rights  reserved, 
dai:  1 0,1 01 6/j.ocemod.2009.02.009 


that  intensification  of  currents  along  a  steep  slope  is  a  manifesta¬ 
tion  of  topographic  Rossby  waves  (TRWs,  Rhines,  1970)  propagat¬ 
ing  along  the  isobaths  such  that  the  shallow  water  is  to  the  right 
(Hamilton,  1990).  The  mechanisms  generating  these  TRWs  are  un¬ 
clear.  The  waves  may  be  excited  by  interactions  of  energetic  meso- 
scale  circulation  features  with  topography.  The  energetic 
circulation  features  are  a  product  of  the  Loop  Current,  a  branch 
of  the  North  Atlantic’s  western  boundary  current  that  dominates 
the  basin-scale  circulation  in  the  Gulf  of  Mexico.  The  Loop  Current 
enters  the  Gulf  through  the  Yucatan  Strait  and  leaves  the  Gulf 
through  the  Straits  of  Florida  (Fig.  la).  The  Loop  Current  flows 
clockwise  with  near-surface  velocities  more  than  0.8  ms1  making 
a  loop  that  extends  northward.  The  Loop  Current  sheds  anticy- 
clonic  mesoscale  eddies  (Loop  Current  eddies)  that  drift  westward 
with  smaller  scale  cyclones,  called  frontal  eddies,  along  their 
periphery  (the  circulation  in  the  Gulf  is  analyzed  in  detail  in  the  lit¬ 
erature,  e  g.,  Sturges  and  Lugo-Fernandez,  2005).  It  is  hypothesized 
that  variations  in  the  position  of  the  front  of  the  Loop  Current,  Loop 
Current  eddy  shedding  events,  and  the  propagation  of  Loop  Current 
eddies  and  their  associated  smaller  cyclones  across  the  Gulf  of 
Mexico  generate  the  TRWs  over  the  deep  continental  slope 
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Fig.  1.  (a)  Bathymetry  of  the  Gulf  of  Mexico  model  domain.  Bathymetry  contours  are  drawn  every  1000  m.  The  Sigsbee  Escarpment  high-resolution  model  domain  is 
indicated  with  a  rectangle,  and  bathymetry  contours  in  this  area  are  drawn  every  500  m.  (b)  Bathymetry  of  the  Sigsbee  Escarpment  model  domain.  Bathymetry  contours  are 
drawn  every  250  m.  The  location  of  the  MMS  12  mooring  is  indicated  in  the  lower  panel 


responsible  for  the  deep  strong  currents  observed  over  the  Sigsbee 
Escarpment  (Hamilton  and  Lugo-Fernandez,  2001;  Oey  and  Lee, 
2002). 

Although  previous  studies  suggest  that  strong  events  over  the 
SE  are  a  manifestation  of  TRWs,  this  idea  has  to  be  revisited  and 
validated  and  alternative  hypotheses  should  be  considered.  For 
example,  recent  observations  suggest  that  strong  currents  along 
the  SE  could  be  related  to  deep  eddies  impinging  upon  the  slope 
(Donohue  et  al.,  2006).  Other  key  questions  concerning  the  deep¬ 
water  dynamics  remain  unanswered,  and  include1  Are  strong  cur¬ 
rents  in  the  deep  Gulf  of  Mexico  dynamically  uncoupled  from  the 
upper-ocean  currents?  What  governs  the  periodicity  of  the  ener¬ 
getic  currents?  Are  wave-like  motions  over  the  deeper  western 


part  of  the  SE  dynamically  related  to  the  upper  SE  further  to  the 
northeast? 

A  suitable  model  of  the  SE  could  help  in  addressing  all  of  these 
questions.  Observational  studies  have  shown  that  the  SE  strongly 
influences  the  local  deep  dynamics  and  resulting  current  patterns 
(Donohue  et  al.,  2006).  Thus,  it  is  envisioned  that  an  appropriate 
model  of  the  SE  should  have  accurate  representation  of  the  bottom 
topography.  The  model  should  have  very  fine  resolution  for  repre¬ 
senting  the  steep  slope  because  the  SE  is  roughly  10-km  wide  at 
the  narrowest  places,  with  vertical  extents  of  500-1000  m.  Due 
to  the  extremely  complex  bottom  topography  in  the  area  and  the 
very  steep  slope  of  the  SE,  developing  a  model  of  this  region  is  chal¬ 
lenging  and  has  numerical  issues  that  must  be  addressed. 
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Several  ocean  models  have  been  used  to  study  the  deep-ocean 
dynamics  in  the  region.  However,  the  models  have  not  had  high 
enough  resolution  as  well  as  sufficiently  accurate  topographic  rep¬ 
resentation  to  simulate  the  near-bottom  strong  currents  and  the 
physical  processes  responsible  for  these  events.  Not  only  has  verti¬ 
cal  resolution  near  the  ocean  bottom  been  too  coarse  in  most  mod¬ 
els,  but  terrain-following  coordinate  models  have  required 
smoothed  topography  and  models  with  geopotential  coordinates 
typically  have  stepped  topography  and  cannot  properly  simulate 
the  propagation  of  the  features  (Dukhovskoy  et  al.,  2006).  Ter¬ 
rain-following  (sigma)  coordinates  permit  accurate  representation 
of  bottom  topography,  but  over  steep  slopes,  sigma  layers  become 
strongly  inclined,  which  causes  numerical  truncation  errors  when 
calculating  horizontal  gradients  in  a  baroclinic  ocean  (Haney, 
1991).  Substantial  efforts  have  been  undertaken  to  overcome  the 
pressure-gradient  problem  in  sigma-coordinate  models.  A  number 
of  approaches  and  numerical  algorithms  that  reduce  this  error 
have  been  developed  including  the  use  of  higher  order  approxima¬ 
tions  (Mihailovic  and  Zanjic,  1984;  Chu  and  Fan,  1997),  weighed 
Jacobian  methods  (Shchepetkin  and  McWilliams.  2003),  reduction 
of  anisotropy  in  sigma-coordinate  models  (Thiem  and  Bemtsen, 
2006),  calculation  of  the  pressure  gradient  from  the  density  inter¬ 
polated  to  z-levels  (Stelling  and  Van  Kester,  1994)  and  others  (see 
overview  in  Ezer  et  al..  2002).  While  some  methods  exhibit  error 
reduction  over  sloping  bottom,  it  is  not  clear  how  they  would  per¬ 
form  over  very  steep  topography. 

In  some  cases  in  order  to  reduce  the  error,  the  bottom  topogra¬ 
phy  is  heavily  filtered  as  in  the  model  studies  by  Oey  (1996)  and 
Oey  and  Lee  (2002)  which  addressed  the  problem  of  TRWs  in  the 
deep  Gulf  of  Mexico.  In  the  latter  study,  the  horizontal  grid  spacing 
of  the  Princeton  Ocean  Model  simulation  used  in  the  experiments 
was  ~5  km  in  the  northern  Gulf  of  Mexico  with  25  sigma  layers. 
The  primary  objective  of  these  studies  was  to  examine  TRWs  over 
gentle  topographic  slopes  in  the  north-central  Gulf  near  the  3000- 
m  isobath.  The  studies  did  not  attempt  to  analyze  currents  over 
steep  regions  such  as  SE  due  to  insufficient  resolution  in  the  model 
(Oey  and  Lee,  2002). 

The  objective  of  this  study  is  to  develop  a  model  of  the  Sigsbee 
Escarpment  capable  of  simulating  dynamical  processes  and  ener¬ 
getic  bottom-intensified  currents  with  characteristics  similar  to 
those  observed.  It  is  envisioned  that  the  methodology  developed 
here  can  eventually  serve  as  a  tool  for  exploring  the  deep-ocean 
dynamics  in  the  region.  For  this  purpose,  the  Navy  Coastal  Ocean 
Model  (NCOM)  with  a  Generalized  Vertical  Coordinate  is  employed 
(Section  3).  A  “vanishing  quasi-sigma”  vertical  coordinate  system 
is  developed  exploiting  the  flexibility  of  the  NCOM  vertical  grid. 
This  vertical  coordinate  system  allows  accurate  representation  of 
topography  in  the  model  with  significantly  reduced  numerical 
truncation  errors  (Section  4).  The  implementation  of  the  new  ver¬ 
tical  grid  is  tested  within  the  SE  domain  prior  to  performing  a  real¬ 
istic  one-year  integration  to  simulate  events  of  high-speed,  near- 
bottom  currents  over  the  escarpment.  The  results  are  analyzed 
and  compared  to  the  observed  characteristics  of  the  energetic 
events  over  the  Sigsbee  Escarpment  (Section  5). 

2.  Background 

2. I.  Description  of  the  region 

The  SE  is  located  at  the  offshore  edge  of  the  Central  Slope  south 
of  the  Louisiana-Texas  Shelf  where  the  seafloor  drops  dramatically 
to  the  Sigsbee  Plain  of  the  Gulf  of  Mexico  Basin  (Fig.  1).  The  local 
maxima  of  the  bottom  slope  in  the  region  define  the  SE,  and  in 
the  model  topography  described  below,  it  locally  exceeds  0.15. 
The  SE  cannot  be  defined  based  on  a  particular  isobath  because 
the  depth  of  the  escarpment  increases  toward  the  southwest.  The 


escarpment  north  of  roughly  26.8°N  and  east  of  90.5°W  drops  from 
around  1400  m  to  a  depth  of  about  2200  m.  To  the  south  of  26.4°N 
and  west  of  approximately  90.8°W,  the  seafloor  drops  abruptly 
from  about  2400  m  to  a  depth  of  3400  m.  For  the  purposes  of  this 
work,  these  depth  ranges  will  define  the  upper  (eastern)  SE  and 
lower  (western)  SE,  respectively  (Fig.  1).  The  SE,  both  upper  and 
lower  regions,  appear  in  maps  of  topographic  contours  as  a  conver¬ 
gence  of  isobaths  to  the  west  of  the  Mississippi  Fan  region. 

2.2.  Historical  observations  of  deep  currents  in  the  Gulf  of  Mexico 

Before  the  1960s,  it  was  generally  assumed  that  very  weak  or  no 
currents  exist  in  the  deep  ocean  (Skoda,  1970).  In  the  late  1950s, 
indirect  evidence  (e.g.,  photographs  of  the  deep-ocean  bottom  with 
ripple  marks,  scour  around  stones,  and  current  lineations  (Hunkins 
et  al.,  1960;  Hollister  and  Heezen,  1966))  suggested  the  possibility 
of  the  existence  of  strong,  near-bottom  currents  in  the  abyss  and 
excited  interest  in  this  topic  among  oceanographers  (e.g.,  Rowe 
and  Menzies,  1968;  Reid,  1969). 

Probably  the  first  evidence  of  strong  bottom  currents  in  the 
deep  Gulf  of  Mexico  was  presented  by  Pequegnat  (1972),  who  re¬ 
ported  substantial  velocities  up  to  0.19  ms1  from  direct  current 
measurements  at  the  Mississippi  Cone  at  depths  of  3000- 
3300  m.  Later,  in  the  mid-1980s,  currents  of  up  to  0.3  m  s  1  were 
measured  over  the  lower  continental  slope  and  deep  basin  of  the 
Gulf  of  Mexico  (Hamilton,  1990).  Hamilton  (1990)  advocates  that 
low-frequency  velocity  fluctuations  observed  in  the  time  series 
could  be  attributed  to  TRWs.  Estimated  periods  of  the  TRWs  are 
25  days  and  40-100  days,  with  wavelengths  of  about  150-250  km. 

Strong  currents  near  the  bottom  over  the  SE  were  found  in 
observations  conducted  by  the  US  Department  of  the  Intenor  Min¬ 
erals  Management  Service  (MMS)  during  a  four-year,  exploratory 
study  of  deepwater  currents  in  the  Gulf  of  Mexico  (Nowlin  et  al., 
2001;  Donohue  et  al.,  2006).  The  observations  were  supplemented 
by  measurements  from  moorings  deployed  in  the  study  area  by  the 
DeepStar  Consortium  and  Louisiana  State  University.  The  depth  of 
the  bottom  current-meters  vaned  from  approximately  10  to  100  m 
above  the  local  bottom.  Two  years  of  current  observations  were 
collected  from  August  1999  to  August  2001  at  the  upper  portion 
of  the  SE  in  the  vicinity  of  Green  Knoll  (Fig.  1b).  These  measure¬ 
ments  are  analyzed  in  Hamilton  and  Lugo-Fernandez  (2001)  and 
Nowlin  et  al.  (2001).  Bottom-intensified  currents  exceeding 
0.85  ms  1  were  recorded  1 1  m  above  the  bottom  at  the  12  mooring 
(27°  13.683'N  89°  58.264'W,  location  is  shown  in  Fig.  lb)  several 
times  during  this  period.  Hamilton  (2007)  indicated  that  the 
upper-  and  lower-level  flows  are  decoupled,  but  speculated  that 
in  some  instances  there  seems  to  be  connection  between  the 
upper-layer  flows  and  the  currents  in  the  lower  layer.  Further 
observations  of  the  SE  region  have  not  provided  conclusive  evi¬ 
dence  of  the  importance  of  a  dynamical  linkage  of  the  upper  and 
lower  layers  over  the  SE  (Donohue  et  al.,  2006). 

During  February  2003-April  2004,  an  array  of  moorings  was  de¬ 
ployed  in  the  area.  Within  this  time  frame,  a  number  of  floats  (neu¬ 
trally  buoyant  and  profiling  floats)  were  launched  in  the  region.  A 
mesoscale-resolving  array  of  27  inverted  echo  sounders  with  pres¬ 
sure  gauges  (Hamilton  et  al.,  2003)  was  deployed  to  enable  a  quan¬ 
titative  mapping  of  the  regional  circulation.  The  observations  have 
revealed  extraordinary  inhomogeneity  in  the  dynamics  over  the  SE. 
Currents  at  the  upper  (eastern)  portion  of  the  SE  differ  in  character 
and  energy  from  currents  at  the  lower  (western)  portion  of  the  SE 
(Donohue  et  al.,  2006).  Near-bottom  mean  currents  (calculated 
from  one-year  records  at  100  m  above  the  bottom)  show  enhanced 
southwestward  flow  along  the  SE.  These  currents  are  persistent  in 
their  direction  and  rarely  reverse.  The  mean  deep  currents  calcu¬ 
lated  from  the  MMS  observations  range  from  0.05  to  0.3  ms  1 
(Fig.  4.2-13,  in  Donohue  et  al.,  2006).  A  particularly  strong  mean 
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current  of  about  0.1  ms'1  is  noticeable  at  the  eastern  part  of  the 
lower  SE  (~26°N  91.5°W).  At  the  upper  SE,  a  local  maximum  of 
the  mean  current  is  observed  at  the  eastern  end  of  the  escarpment 
where  the  isobaths  converge  (27-27.3°N,  89.2-89.9°W).  Over  the 
shallower  depths  to  the  north  of  the  escarpment,  the  mean  cur¬ 
rents  are  mostly  weaker  revealing  the  insulating  effect  of  the  steep 
slope  (Donohue  et  al.t  2006).  Estimated  from  the  observations,  the 
standard  deviation  (STD)  ellipses  of  the  mean  current  vectors  along 
the  SE  have  major  axes  that  are  rotated  at  small  angles  to  the  iso¬ 
baths.  Over  the  less  steep  slope,  the  ellipses  are  more  circular  and 
the  major  axes  are  oriented  at  greater  angles  to  the  isobaths. 

Extreme  events  over  the  escarpment  are  episodic  in  time  and 
localized  in  space.  Analysis  of  the  measurements  suggests  that  dur¬ 
ing  2003-2004  the  deep  currents  were  in  the  upper  quartile  of 
their  speed  range  in  only  0.5- 1.5%  of  the  records  (Donohue  et  al., 
2006).  Regions  of  maximum  and  minimum  speed  of  the  currents 
can  be  close  to  each  other  along  the  escarpment.  Hamilton 
(2007)  has  noticed  drastically  different  current  speeds  measured 
at  the  upper  SE  at  two  moorings  6  km  apart.  At  the  12  mooring 
location  (shown  in  Fig.  1b),  recorded  currents  at  11  m  above  the 
bottom  exceeded  0.7  m  s  \  with  speeds  occasionally  approaching 
1  ms1.  The  speeds  at  the  other  mooring  located  6  km  to  the  west 
along  the  slope  were  approximately  half  of  those  at  the  12  location. 
During  March  2003-March  2004,  the  maximum  speeds  observed 
in  the  area  were  not  as  extreme  as  those  measured  by  the  Septem¬ 
ber  1999  12  mooring.  The  maximum  currents  measured  at  instru¬ 
ments  from  100  to  500  m  above  the  bottom  were  mostly 
concentrated  in  two  locations:  one  at  the  upper  SE  (~27-27.5°N, 
89.2-90.2°W)  with  speeds  of  0.5-0.6  ms  \  and  the  other  at  the 
lower  SE  (~26-26.3°N,  90.8-9 1.3°W)  with  currents  of  0.4- 
0.5  ms  \  Intensification  of  the  currents  in  the  lower  part  of  the 
water  column  as  reported  in  Hamilton  and  Lugo-Fernandez 
(2001 )  has  also  been  observed  during  these  observations  (Hamil¬ 
ton,  2007).  In  these  current  profiles,  however,  a  reduction  of  speed 
near  the  bottom  is  evident  and  likely  due  to  a  frictional  bottom 
boundary  layer  resulting  in  the  maximum  speeds  occurring  at 
some  distance  above  the  bottom.  For  example,  the  greatest  mean 
speed  in  the  deep  layer  measured  at  a  mooring  over  the  lower  SE 
(~91°W,  26°N)  was  200  m  above  the  bottom  (Donohue  et  al , 
2006). 

Spectral  analyses  of  the  velocity  time  series  have  identified  the 
ranges  of  predominant  periodicities  of  the  low-frequency  oscilla¬ 
tions  over  the  SE.  Kinetic  energy  spectra  of  the  velocity  records 
do  not  have  a  well  defined  peak,  but  the  energy  is  concentrated 
within  a  wide  frequency  range.  Periodicities  of  the  fluctuations  de¬ 
pend  on  the  location  of  the  records,  both  the  east-west  position 
along  the  escarpment  and  the  depth.  Typically,  the  estimated  time 
scales  of  the  deep,  low-frequency  motions  in  the  region  vary  from 
10  to  40  days  (Hamilton  and  Lugo-Fernandez,  2001;  Donohue 
et  al.,  2006;  Hamilton,  2007). 

2.3.  Prominent  features  of  the  deep-ocean  dynamics  over  the  SE 

Based  on  the  historical  observations,  a  fundamental  set  of 
prominent  features  are  identified.  It  is  expected  that  for  a  model 
to  be  deemed  useful  for  studying  these  dynamics,  it  should  be  able 
to  reproduce  these  characteristics  of  the  deep  currents  in  the  re¬ 
gion  about  the  SE: 

Mean  currents.  Along  the  SE  the  near-bottom  mean  currents  are 
directed  dominantly  southwestward.  Speeds  of  the  mean  currents 
range  from  0.05  to  0.3  m  s'1.  The  mean  currents  are  weaker  over 
the  shallower  (upper)  part  of  the  SE.  Along  the  slope,  intensifica¬ 
tion  of  the  mean  currents  is  observed  at  the  eastern  sides  of  the 
lower  and  upper  SE.  The  major  axes  of  the  STD  ellipses  are  oriented 
at  small  angles  to  the  isobaths  over  the  escarpment. 


Extreme  events.  Episodic  extreme  events  affect  the  SE  with 
maxima  occurring  just  below  or  on  the  escarpment  in  two  separate 
regions,  one  near  the  eastern  end  of  the  upper  SE  (~89.5-90°W, 
27-27.3°N)  and  the  other  over  the  eastern  lower  SE  (91-91.2°W, 
26-26.2°N).  At  the  upper  SE,  speeds  of  the  observed  extreme  cur¬ 
rents  have  been  measured  in  the  range  of  0.6  to  greater  than 
1ms  \  At  the  lower  SE,  observed  maximum  speeds  over 
0.65  m  s  1  have  been  observed. 

Time  scales  of  deep  currents.  The  estimated  periodicity  of  the 
low-frequency  oscillations  in  the  velocity  records  over  the  SE  is 
10-40  days  with  some  longer  periods  up  to  100  days. 


3.  Numerical  model  description  and  simulations 

3. 1.  The  Navy  Coastal  Ocean  Model 

All  numerical  simulations  used  for  this  project  are  applications 
of  the  Navy  Coastal  Ocean  Model  (NCOM)  (Martin,  2000),  devel¬ 
oped  at  the  U.S.  Naval  Research  Laboratory  for  operational  applica¬ 
tions.  Specifically,  the  NCOM  Version  3.2  code  is  used.  The  NCOM  is 
a  primitive-equation,  three-dimensional,  ocean  model  with  the 
hydrostatic  and  Boussinesq  approximations.  The  model  uses  an 
Arakawa  C  staggered  grid.  It  has  several  numerical  differencing 
and  integration  options.  For  this  study,  second-order-centered 
averages  and  differencing  methods  are  used  for  spatial  interpola¬ 
tion  and  approximation  of  spatial  derivatives,  with  the  leapfrog 
semi-implicit  method  used  for  time  integration.  The  bottom  fric¬ 
tion  used  in  the  simulation  is  parameterized  with  a  quadratic  drag 
law  and  implicitly  calculated  by  the  model.  The  drag  coefficient  cb 
is  calculated  as 

cb  -  max  <  777  t  ,cbmin  1  (1) 

M&)  1 

Where,  K  is  Von  Karman’s  constant,  cbmjn  is  the  minimum  value 
that  the  bottom  drag  coefficient  can  attain  (set  at  0.0025  in  this 
study),  A zb  is  the  height  of  the  bottom  grid  cell  and  z0  the  bottom 
roughness  (in  this  study  z0  «  0.01  m). 

The  NCOM  belongs  to  a  class  of  hybrid-vertical-coordinate 
ocean  models.  Early  versions  of  the  NCOM  combine  geopotential 
(z-level)  coordinates  with  terrain-following  (sigma)  coordinates. 
The  sigma  coordinates  are  specified  as  a  fraction  of  the  total 
water-column  depth  from  the  free  surface  to  either  the  ocean  bot¬ 
tom,  or  in  the  case  of  the  NCOM’s  hybrid  vertical  grid,  some  spec¬ 
ified  depth  level  if  the  ocean  bottom  is  deeper.  The  a  coordinates 
need  not  be  equally  spaced  and  can  be  stretched.  The  z-levels  are 
at  fixed  depths  below  a  resting  ocean  surface  (geopotentia!)  and 
the  ocean  model  topography  is  rounded  at  each  horizontal  point 
to  the  nearest  z-level.  The  use  of  z-levels,  while  eliminating  the 
pressure-gradient  error  of  sigma-coordinate  models  (Mellor  et  al., 
1994,  1998),  results  in  a  step-like  representation  of  topography 
and  presents  complications  when  higher  resolution  is  desired  near 
the  ocean  bottom.  Dukhovskoy  et  al.  (2006)  showed  that  topo¬ 
graphic  Rossby  waves  cannot  be  simulated  properly  in  a  z-level 
model  unless  a  large  number  of  levels  are  used  so  that  the  size  of 
the  topographic  steps  are  reduced  to  form  a  reasonably  accurate 
representation  of  the  topography. 

NCOM  Versions  3.1  and  later  feature  a  new  extension  to  the  hy¬ 
brid  vertical  grid,  termed  by  the  model  developers  as  a  Generalized 
Vertical  Coordinate  (GVC).  This  new  vertical  coordinate  system  of¬ 
fers  traditional  sigma  coordinates  (that  can  vary  in  time  as  the  total 
water-column  depth  changes  due  to  free-surface  deviations),  fixed 
depth  (or  static)  coordinates  with  vertical  distributions  varying  as 
a  function  of  horizontal  location,  traditional  z-level  coordinates,  or 
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z-level  grids  with  shaving  of  the  bottom  grid  cell  to  match  the  real¬ 
istic  topography.  Additionally,  grid  cells,  even  on  the  sigma  grid, 
can  be  masked  out  so  they  are  not  active  grid  cells.  The  generalized 
nature  of  the  vertical  coordinate  system  allows  for  great  flexibility, 
which  has  been  exploited  for  this  study  through  the  development 
of  what  will  be  termed  vanishing  quasi-sigma  (VQS)  coordinates. 

3.2.  Vanishing  quasi -sigma  vertical  coordinate 

The  VQS  grid  is  defined  as  a  function  of  the  local  depth  of  the 
ocean  bottom.  As  the  ocean  becomes  shallower,  the  number  of  ac¬ 
tive,  sigma-type  coordinates  is  reduced.  That  is,  the  deepest  coor¬ 
dinates  "vanish”  as  they  are  masked  out  (Fig.  2a  and  d).  The  grid 
is  generated  as  follows.  A  sequence  of  reference  vertical  grids  is 
generated  for  specified  (increasing)  bottom  depths.  Each  reference 
vertical  grid  is  computed  using  a  uniform  (linear)  stretching  func¬ 
tion  to  provide  higher  resolution  near  the  bottom  and  the  surface 
(Fig.  2b).  The  stretching  parameter  is  defined  through  interval 
halving  iteration  until  the  thickness  of  the  bottom  and  surface  lay¬ 
ers  do  not  exceed  desired  values  (<15  m  in  this  study).  The  refer¬ 
ence  vertical  grids  have  different  numbers  of  vertical  layers 
(increasing  with  depth)  for  given  depths  and  specified  stretching 
factors.  For  a  model  grid  point  (l,  j )  with  a  specific  depth  /i(i,  j) 
the  two  reference  grids  corresponding  to  reference  depths  above 


and  below  h(\,j)  are  identified.  The  vertical  grid  is  computed  for 
that  point  using  linear  weighing  of  the  two  reference  grids.  The 
minimum  bottom  layer  vanishes  when  it  becomes  thinner  than  a 
specified  threshold  value  (3  m  in  this  study).  The  generated  grid 
is  forced  to  maintain  a  minimum  number  of  layers  (10  in  this 
study)  in  shallow  regions,  where  the  layer  thicknesses  may  be  less 
than  3  m. 

The  VQS  scheme  has  the  benefit  of  allowing  for  the  vertical  grid 
to  closely  follow  the  ocean  bottom  slope,  but  restricts  the  depth 
range  for  a  given  coordinate  surface  reducing  its  slope  near  steep 
topographic  features  when  compared  to  the  classic  sigma-coordi¬ 
nate  application  (Fig.  2a  and  c).  Reducing  the  slopes  of  the  sig- 
ma-coordinate  surfaces  reduces  the  numerical  truncation  errors 
(demonstrated  in  the  following  sections),  which  in  classic  sigma 
models  can  be  done  by  reducing  the  slope  of  the  escarpment  (with¬ 
out  increasing  the  horizontal  resolution). 

3.3.  Application  of  the  NCOM  to  the  Sigsbee  Escarpment  region 

To  achieve  the  very  high  resolution  necessary  for  properly 
resolving  the  SE  and  simulating  the  impacts  of  larger-scale  circula¬ 
tion  features  on  the  region,  a  two-level  nested  modeling  system  is 
applied  using  the  NCOM.  Two  model  domains  are  configured  for 
this  nested  approach  using  the  VQS  vertical  grid.  The  outer  domain 


Fig.  2.  (a)  Vertical  coordinate  surfaces  along  90  W  from  the  high-resolution  SE  model  using  the  VQS  coordinate  system.  Every  10th  coordinate  surface  is  highlighted  for 
clarity.  A  rectangle  indicates  the  region  shown  in  (d).  (b)  Vertical  grid  spacing  (m)  versus  depth  at  three  locations  (shown  by  bullets  in  (a)),  (c)  An  evenly  stretched  traditional 
sigma  grid  with  80  layers  is  shown  (used  in  test  runs).  Note  in  the  traditional  grid  strong  inclination  of  sigma  surfaces  over  steep  slopes  through  the  water  column  as  well  as 
overburdening  of  the  shallow  region  with  sigma  layers,  (d)  Magnification  of  the  VQS  coordinate  grid  bounded  by  the  rectangle  in  (a). 
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is  a  1/20°  horizontal  resolution  (in  latitude  and  longitude  between 
like  variables  on  the  Arakawa  C  grid)  encompassing  the  entire  Gulf 
of  Mexico  and  northwestern  Caribbean  (GoM  domain.  Fig.  la). 
Next  is  a  SE  domain  with  a  30  arc  second  (1/120°)  grid  (Fig.  lb) 
and  is  focused  on  the  region  surrounding  the  escarpment.  The  out¬ 
er  GoM  model  domain  is  integrated  for  several  (10)  years  to  simu¬ 
late  the  Loop  Current  and  the  mesoscale  eddy  field  of  the  GoM  and 
to  choose  an  appropriate  case  for  a  high-resolution  nested  model 
simulation  of  the  SE.  One  particular  year  with  intensified  Loop  Cur¬ 
rent  Eddy  activity  over  the  SE  region  in  the  GoM  model  domain  is 
identified  as  a  case  for  which  to  run  the  nested  SE  model 
experiment. 

3.3.1.  Gulf  of  Mexico  domain 

Circulation  of  the  Gulf  of  Mexico  and  northwestern  Caribbean  is 
simulated  in  a  model  domain  encompassing  the  region  98.15°W  to 
80.6°W  and  15.55°N  to  31.5°N  (Fig.  la).  The  eastern  boundary  is 
open  in  the  Caribbean  Sea  and  Florida  Strait,  and  a  small  part  of 
the  southern  boundary  extending  over  the  shallow  Nicaragua  Shelf 
has  been  artificially  closed.  The  northern  and  western  boundaries 
are  naturally  closed.  The  topography  for  the  model  is  derived  from 
the  NOAA  ET0P02  2-min-resolution  data  set  with  some  regional 
editing  and  interpolated  to  the  model  gnd  before  applying  a  9- 
point  box  average  twice.  The  simulation  domain  and  forcing  are 
based  on  that  described  in  Morey  et  al.  (2003). 

The  vertical  grid  is  configured  as  follows.  Below  the  free  surface 
(the  first  vertical  grid  level),  10  static  sigma  layers  are  evenly 
spaced  in  the  upper  150  m  Thus  the  top-most  layer  varies  its 
thickness  with  the  free  ocean  surface.  Below  1 50  m,  the  VQS  grid 
is  applied  with  a  maximum  of  50  additional  layers  (yielding  a  total 
of  up  to  60  vertical  layers). 

The  model  is  initialized  at  rest  with  temperature  and  salinity 
interpolated  from  the  World  Ocean  Atlas  1998  (WOA98)January  cli¬ 
matology  (Conkright  et  al.,  1998).  Surface  boundary  conditions  for 
momentum  and  temperature  are  obtained  from  the  DaSilva  et  al. 
(1994)  monthly-climatological  heat  and  radiative  fluxes  and  wind 
stresses.  A  salinity  flux  is  prescribed  to  be  constant  in  space  and  time 
to  balance  the  annual  average  fresh  water  input  by  the  thirty  rivers 
in  the  model  domain.  The  rivers  are  treated  as  mass  sources  of  fresh 
water  with  flow  rates  derived  from  monthly-climatological  USGS 
data  or  other  sources  for  Mexican  rivers.  Open  boundary  conditions 
are  obtained  by  adjusting  the  baroclinic  flow  field  at  the  boundary  to 
the  WOA98  average  temperature  and  salinity  field.  This  forces  the 
incoming  Caribbean  Current  and  exiting  Florida  Current,  producing 
a  mean  Yucatan  Channel  transport  of  26  Sv. 

Integrating  the  model  from  its  initial  state,  a  Loop  Current  is 
rapidly  formed  The  first  Loop  Current  Eddy  is  shed  in  the  sixth 
month  of  integration  The  model  is  integrated  for  10  years  for  anal¬ 
ysis  The  Loop  Current  penetration  and  eddy  shedding  behavior  is 
similar  to  observations  reported  in  prior  studies.  For  example, 
the  average  period  between  Loop  Current  Eddy  separations  from 
8  years  of  model  data  is  10.2  months,  varying  between  4  and 
22  months.  This  can  be  compared  to  the  Sturges  and  Leben 
(2000)  estimate  from  satellite  observations  of  9.5  months,  varying 
between  3  and  17  months,  and  Leben’s  (2005)  estimate  of 
9.4  months  from  16  altimetry-observed,  eddy-separation  events. 

3.3.2.  High -resolution  Sigsbee  Escarpment  domain 

A  second  model  domain  is  configured  to  include  the  SE  region 
and  surrounding  Gulf  of  Mexico  waters  from  95°W  to  88°  30.5'W 
and  25°N  to  28°29.5'N  (Fig.  lb).  The  horizontal  resolution  of  the 
model  grid  is  30  arc  seconds  (1/120°).  Topography  is  derived  from 
the  NOAA  NODC  30-arc-second  gridded  Coastal  Relief  and  Digital 
Elevation  Map,  with  a  nine-point  centered-average  filter  applied 
three  times  to  the  bathymetric  data  set  prior  to  use  in  the  model. 
The  bathymetry  is  blended  with  the  coarser-resolution  GoM  1/ 


20°  model  topography  near  the  edges  of  the  domain  so  that  the 
two  topographies  match  at  the  boundaries.  The  vertical  grid  has 
a  maximum  of  80  layers  and  is  configured  with  the  VQS  grid  (de¬ 
scribed  in  Section  3.2)  (Fig.  2a).  The  thickness  of  the  bottommost 
grid  cell  is  between  3  and  1 5  m 

The  high-resolution  SE  model  is  nested  within  the  larger-scale 
GoM  model  and  run  for  one  year  selected  from  the  10-year  GoM 
integration.  Initial  fields  (temperature,  salinity,  currents,  and  sea- 
surface  elevation)  and  boundary  conditions  are  interpolated  from 
the  GoM  model  to  the  fine-resolution  gnd.  Coupling  is  one-way 
in  that  outgoing  information  is  simply  radiated  away  and  does 
not  feed  back  into  the  large-scale  model.  For  this  realistic  simula¬ 
tion,  the  same  surface  fluxes  applied  to  the  GoM  model  domain 
are  used  to  force  the  nested  SE  model  domain.  Prior  to  conducting 
the  realistic  simulation  model  experiments,  test  model  runs  are 
performed  with  the  SE  model  domain  as  described  in  Section  4.2 
to  evaluate  the  performance  of  the  new  VQS  grid. 


4.  Evaluation  of  the  VQS  vertical  coordinate  system 


4. 1.  Sigma-coordinate  pressure-gradient  errors 


The  manner  in  which  the  new  vertical  grid  is  defined  distin¬ 
guishes  it  from  the  classic  sigma-coordinate  system  while  preserv¬ 
ing  certain  useful  features  of  the  sigma  coordinate.  The 
conventional  definition  of  sigma  (<r)  coordinates  in  two  dimensions 
is 


z  -  n(*  0 

H(x)itj(x,t) 


(2) 


where  (x,  z)  are  the  horizontal  and  vertical  Cartesian  coordinates,  H 
is  the  water-column  depth  below  a  resting  free  surface  (z  -  0),  and  t] 
is  the  free-surface  deviation  at  time  t.  Coordinate  transformation  (2) 
defines  surfaces  such  that  the  bottom  surface  follows  the  terrain 
(thus  called  “terrain-following”)  and  the  upper  surface  matches 
the  sea  surface.  In  a  sigma-coordinate  ocean  model,  sigma  surfaces 
define  the  upper  and  lower  interfaces  of  volume  elements  within 
the  sigma  layers,  but  their  depths  can  change  in  time  with  varia¬ 
tions  in  the  sea  surface.  Hence,  the  vertical  model  layers  coincide 
with  the  vertical  sigma  layers  (Fig.  2c).  Within  a  fixed  vertical  layer 
k,  the  sigma  coordinate  does  not  change  horizontally: 


(3) 


In  the  VQS  coordinate  system,  as  in  classic  sigma  coordinates, 
the  bottom  surface  follows  the  terrain,  the  top-most  surface  fol¬ 
lows  the  free  surface  of  the  ocean,  and  the  depths  of  interior  coor¬ 
dinate  surfaces  are  computed  as  fractions  of  the  total  water- 
column  depth.  These  fractions,  however,  are  not  fixed  over  the  do¬ 
main  but  are  allowed  to  vary  at  each  horizontal  gnd  point,  and  the 
total  number  of  vertical  levels  also  varies  with  the  seafloor  depth. 
Thus,  the  VQS  coordinate  is  defined  algorithmically  as  described  in 
Section  3.2  as  opposed  to  a  simple  analytical  expression  like  (2).  A 
key  difference  is  that  in  the  VQS  vertical  coordinate  system,  the 
vertical  coordinate  surfaces  can  intersect  surfaces  of  constant  sig¬ 
ma  coordinates  and  (3)  no  longer  applies  (compare  Fig.  2a  and  c). 

Fora  two-dimensional  case,  the  horizontal  gradient  of  pressure 
in  Cartesian  coordinates,  when  evaluated  on  the  VQS  grid  (which 
are  not  level  surfaces),  is 

ap\  =  <^\  _dpM\ 
dx)  z  dx)  k  dkdxj  z 

On  VQS  coordinates,  as  in  sigma  coordinates,  the  horizontal  gra¬ 
dient  of  pressure  is  evaluated  as  the  pressure  gradient  along  a  fixed 
kth  layer  (the  first  term  on  the  right-hand  side  of  (4))  and  the 
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“hydrostatic  correction”  (the  second  term),  which  is  proportional 
to  the  slope  of  the  vertical  layers.  Analogous  to  the  computation 
on  sigma  coordinates,  the  truncation  error  in  the  calculation  of 
the  pressure  gradient  is  a  function  of  the  slope  of  the  layers.  When 
the  layers  are  horizontal,  the  second  term  is  zero  and  so  is  the  er¬ 
ror.  Thus,  the  smaller  slopes  of  the  vertical  layers  in  the  VQS  grid 
should  lead  to  smaller  truncation  errors  in  the  calculation  of  the 
horizontal  pressure  gradient. 

4.2.  Model  test  runs 

Several  metrics  have  been  suggested  to  estimate  truncation  er¬ 
rors  in  sigma  models  (Haney,  1991;  Mellor  et  al.,  1994,  1998).  The 
utility  of  Haney’s  parameter  (1991)  has  been  critiqued  in  Mellor 
et  al.  (1994, 1998)  that  instead  suggest  analytically  derived  metrics 
for  Sigma  Errors  of  the  First  and  Second  Kinds  (SEFK  and  SESK, 
respectively).  The  analytical  forms  of  these  metrics  are  not  applica¬ 
ble  to  the  VQS  grid  because  this  grid  is  not  truly  based  on  sigma 
coordinates  and  an  analytical  relation  between  sigma  surfaces 
and  vertical  layers  is  not  tractable.  To  evaluate  the  performance 
of  VQS  coordinates  within  a  numerical  model  and  to  compare  with 
a  classic  sigma  grid,  idealized  numerical  tests  using  a  suite  of  test 
model  runs  within  the  high-resolution  SE  domain  described  in  Sec¬ 
tion  3.3  are  conducted. 

For  the  test  experiments,  the  NCOM  SE  model  domain  is  config¬ 
ured  with  three  different  vertical  grids:  80  evenly  spaced  non-van¬ 
ishing  (or  classic)  sigma  (NVS)  layers  computed  by  (2),  40  evenly 
spaced  NVS  layers,  and  80  VQS  layers  computed  as  described  in 
Section  3  2.  Hence,  these  experiments  allow  comparison  of  the  per¬ 
formance  of  the  new  VQS  vertical  grid,  which  is  used  for  realistic 
simulation  of  the  deep  currents  in  the  following  sections,  to  the 
classic  sigma  vertical  grid.  Each  model  is  initialized  at  rest  with 
horizontally  homogeneous  density  fields  derived  from  long-term- 
average  density  profiles  for  the  region  from  Conkright  et  al 
(1998).  The  experiments  are  run  prognostically  for  20  days  with 
no  forcing  at  the  surface  or  boundaries.  The  theoretical  solution 
is  that  the  ocean  remains  at  rest  with  no  change  in  the  density 
field. 


4.2.1.  Slopes  of  vertical  coordinate  surfaces 

In  general,  the  slopes  of  the  near-bottom  vertical  layers  in  the 
model  with  VQS  coordinate  are  noticeably  less  steep  over  the  SE 
than  in  the  model  with  80  and  40  NVS  layers  which  has  similar 
slopes  of  the  sigma  layers  (Fig.  3,  40  NVS  case  is  not  shown).  Over 
the  upper  SE.  the  VQS  grid  has  remarkably  reduced  slopes  of  the 
vertical  layers  with  a  maximum  slope  of  0.1  (Fig.  3d),  which  is 
one-third  smaller  than  in  the  NVS  grids  (Fig.  3b),  but  the  maximum 
slopes  in  the  VQS  grid  (Fig.  3c)  are  nearly  the  same  as  in  the  NVS 
grid  over  the  lower  SE  (Fig.  3a).  In  the  VQS  model,  however,  large 
slopes  are  confined  to  several  grid  cells  near  the  bottom  over  the 
steepest  part  of  the  escarpment,  while  in  the  NVS  models,  large 
slopes  are  found  throughout  the  water  column  over  a  wider  range 
of  the  sloping  bottom.  Thus,  in  the  VQS  grid,  the  upper  layers  are 
not  substantially  inclined  due  to  the  steep  bathymetry  underneath 
and  the  coordinate  surfaces  are  nearly  horizontal  in  the  upper 
ocean  where  the  strong  vertical  density  gradients  occur  (compare 
Fig.  2a-c). 

As  discussed  in  a  number  of  papers  (e.g.,  Auclair  et  al.,  2000), 
sharp  density  gradients,  when  associated  with  inclined  vertical 
coordinate  surfaces,  lead  to  large  truncation  errors  in  the  computa¬ 
tion  of  the  pressure-gradient  term.  In  the  NVS  models,  sigma  layers 
preserve  distinguishably  non-zero  slopes  to  above  1 000  m  over  the 
upper  SE.  Doubling  the  number  of  vertical  layers  in  NVS  models 
does  not  reduce  the  slopes  of  the  vertical  layers,  suggesting  that 
addition  of  sigma  layers  does  not  reduce  the  error,  in  agreement 
with  Mellor  et  al.  (1994). 

Parameters  in  the  algorithm  generating  the  VQS  grid  are  chosen 
such  that  the  layers  vanish  more  rapidly,  providing  smaller  slopes, 
over  the  upper  SE  than  over  the  lower  SE.  This  choice  is  dictated  by 
the  assumption  that  regions  with  weaker  stratification  (lower  SE) 
can  better  tolerate  greater  slopes  of  the  vertical  layers  than  regions 
with  stronger  stratification  (upper  SE)  (Auclair  et  al.,  2000). 

4.2.2.  Density  anomalies 

From  (4),  truncation  errors  depend  on  the  pressure  gradient 
across  the  vertical  coordinate  surfaces  dp/dk  (which,  in  the  case 
of  classic  sigma  coordinates,  is  equivalent  to  dp /do).  At  the  resting 
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Fig.  3.  Cross-sections  of  slopes  of  vertical  coordinate  surfaces.  The  cross-sections  are  at  the  lower  (left  column)  and  upper  (right  column)  SE  along  90.97°W  and  89.93°W, 
respectively.  The  upper  row  (a  and  b)  is  from  the  model  with  80  non-vanishing  sigma  layers,  and  the  bottom  row  (c  and  d)  is  for  80  VQS  layers  (configuration  used  for 
simulating  the  deep  currents  over  SE). 
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initial  state,  the  LHS  of  (4)  is  identically  zero,  so  the  first  term  on 
the  RHS  of  (4)  should  be  balanced  by  the  second  term  on  the 
RHS.  Due  to  the  truncation  error,  however,  this  balance  is  violated 
as  the  model  integrates.  After  the  first  time  step,  the  horizontal 
pressure  gradient  is  no  longer  zero  and  the  initially  motionless 
ocean  starts  moving.  Due  to  the  resulting  mixing  and  advection, 
the  model  density  field  becomes  no  longer  horizontally 
homogeneous. 

Cross-sections  of  the  density  anomaly  (the  initial  density  field  is 
subtracted)  after  20  days  of  integration  (Fig.  4,  the  40  NVS  case  is 
not  shown  since  it  is  similar  to  the  80  NVS  result)  show  that  the 
density  field  in  all  the  models  has  drifted  from  the  initial  density 
distribution  and  noticeable  horizontal  gradients  have  developed 
throughout  the  water  column.  In  the  models  with  80  (Fig.  4a) 
and  40  NVS  layers,  strong  horizontal  density  gradients  develop 
over  the  lower  SE  to  a  height  of  1000-1500  m  above  the  bottom. 
The  experiment  with  the  VQS  grid  (Fig.  4c)  preserves  better  hori¬ 
zontal  homogeneity  of  the  density  fields,  although  weak  anomalies 
are  developed  near  the  lower  SE  slope.  At  the  upper  SE,  the  density 
fields  of  the  NVS-grid-model  experiments  reveal  strong  horizontal 
inhomogeneity  with  maximum  anomalies  developed  over  the 
escarpment  (Fig.  4b),  The  model  with  VQS  layers  has  weaker  hori¬ 
zontal  density  gradients  over  the  upper  SE  with  remarkably  small 
gradients  at  its  steepest  part  (Fig.  4d).  It  is  noteworthy  that  in  the 
model  with  the  NVS  grid,  the  maximum  density  anomalies  occur  in 
the  surface  layer,  where  the  VQS  grid  model  develops  only  small 
anomalies  during  the  test  run. 

4.2.3.  Maximum  speeds 

In  the  models  with  the  NVS  grid,  non-zero  velocities  arising 
from  truncation  errors  in  the  horizontal  pressure-gradient  calcula¬ 
tion  are  largest  near  the  upper  and  lower  SE  (Fig.  5a  and  b,  only  80 
NVS  case  is  shown).  After  20  days  of  the  simulation,  the  maximum 
speed  over  the  lower  SE  (0  08-0.10  m  s  1 )  is  roughly  one-third  of 
the  maximum  speed  over  the  upper  SE  (0.25  ms  1 )  in  these  exper¬ 
iments.  This  is  because  the  stratification  at  the  upper  SE  is  stronger, 
resulting  in  larger  truncation  errors  and,  thus,  larger  density  anom¬ 
alies,  while  the  slopes  at  the  cross-sections  shown  are  comparable. 


Over  the  upper  SE  in  the  NVS-gnd  experiments,  these  spurious 
velocities  have  magnitudes  roughly  2.5  time  greater  than  the  mean 
near-bottom  speed  measured  in  this  area  (~0. 1  ms1)  and  reach 
nearly  one-half  of  the  maximum  deep  current  (~0.5-0.6  ms  1 )  ob¬ 
served  at  the  MMS  moorings  over  the  upper  SE  in  2003-2004 
(Donohue  et  al„  2006),  At  the  lower  SE,  speeds  resulting  from  trun¬ 
cation  errors  (0.08-0.1  ms  *)  are  comparable  to  the  observed 
mean  (~0.1  ms  l).  Thus,  it  would  be  problematic  to  use  this  model 
with  sigma  coordinates  to  produce  accurate  simulations  of  the 
deep  currents  over  the  SE. 

Large  velocities  in  the  NVS-grid-model  experiments  appear  in 
the  lower  water  column  from  the  upper  SE  through  the  eastern 
part  of  the  lower  SE  (Fig.  6a  and  b).  It  is  worth  mentioning  that 
increasing  the  number  of  vertical  layers  does  not  reduce  the  erro¬ 
neous  velocities,  as  the  maps  of  speeds  near  the  bottom  look  sim¬ 
ilar  in  the  models  with  80  and  40  sigma  layers.  These  velocities  are 
a  manifestation  of  the  case  when  the  Sigma  Error  of  the  First  Kind 
(Mellor  et  al.,  1994)  is  dominated  by  the  change  in  water-column 
depth  between  grid  points  and  is  not  sensitive  to  the  number  of 
vertical  layers. 

The  smaller  slopes  of  the  vertical  coordinate  surfaces  in  the  VQS 
grid  result  in  much  weaker  and  more  localized  deep  currents  over 
the  lower  and  upper  SE  (Figs.  5c  and  d,  6c)  compared  to  the  NVS 
cases.  Over  the  upper  SE,  speeds  are  less  than  0.1  m  s  \  except 
for  at  several  grid  points  very  close  to  the  escarpment  where 
speeds  reach  a  maximum  of  0.11-0.12  ms  In  the  VQS  model, 
spurious  velocities  remain  very  small  in  amplitude  (<0.02  ms  *) 
across  most  of  the  slope  of  the  lower  SE,  with  a  maximum  speed 
of  0.05  m  s  1  occurring  600  m  above  the  shallower  part  of  the  low¬ 
er  SE. 

Time  series  of  the  maximum  (Fig.  7,  upper  row)  and  spatially- 
averaged  (Fig.  7,  bottom  row)  speeds  are  computed  over  the  east¬ 
ern  parts  of  the  lower  (91.5W-90.8W)  and  upper  (90.5W-89.8W) 
SE  for  the  deepest  layer  at  each  grid  point  for  the  three  model 
experiments.  Domains  for  spatial  averaging  are  defined  by  the 
1500  and  2000  m  isobaths  at  the  upper  SE,  and  the  2500 
and  3000  m  isobaths  along  the  lower  SE.  These  metrics  show 
clearly  the  utility  of  the  VQS  grid.  The  maximum  speeds  resulting 
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Fig.  4.  Cross-sections  of  density  anomalies  developed  from  day  0  to  day  20  of  the  test  run  in  the  models  with  different  configuration  of  the  vertical  coordinates.  The  cross- 
sections  are  at  the  lower  (left  column)  and  upper  (right  column)  SE  along  90.97W  and  89.93.  respectively  The  upper  row  (a  and  b)  is  the  model  with  80  NVS  layers,  and  the 
bottom  row  (c  and  d)  is  from  the  model  with  80  VQS  layers. 
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Fig-  5.  Cross-sections  of  speeds  resulting  from  truncation  errors  in  the  calculation  of  horizontal  pressure  gradient  in  the  models  with  different  configuration  of  the  vertical 
coordinates.  The  cross-sections  are  at  the  lower  (left  column)  and  upper  (right  column)  SE  along  90.97Wand  89.93W.  respectively.  The  upper  row  (a  and  b)  is  the  model  with 
80  NVS  layers,  and  the  bottom  row  (c  and  d)  is  from  the  model  with  80  VQS  layers.  Note  different  scales  for  the  Lower  and  Upper  SE.  ’‘Max  |U|"  indicates  values  of  maximum 
speed  over  the  Escarpment. 


from  truncation  errors  are  several  times  smaller  in  the  model  with 
the  VQS  gnd  than  in  the  NVS-grid  models.  The  difference  is  more 
noticeable  over  the  upper  SE  (Fig.  7b)  where  the  maximum  speed 
in  the  model  with  the  VQS  is  about  0.05  m  s  1  while  in  the  models 
with  NVS  grids  the  maximum  speed  is  0.15-0.18  ms  1  for  80  lay¬ 
ers  and  0.21-0.23  ms  1  for  40  layers. 

An  interesting  detail  seen  in  Fig.  7  is  slower  spurious  currents 
simulated  in  the  bottommost  grid  cell  of  the  80-layer  NVS-grid 
model  compared  to  the  40-layer  model.  This  difference  is  not 
noticeable  in  the  plot  of  the  maximum  speed  computed  over  the 
200-m  thick  near-bottom  layer  (Fig.  6).  A  possible  explanation  is 
that  the  bottommost  grid  cell  in  the  80-layer  model  is  thinner 
and  does  not  extend  as  high  off  the  bottom  as  the  bottommost  grid 
cell  in  the  40-layer  model.  Note  that  the  drag  coefficient  is  inver¬ 
sely  proportional  to  the  logarithm  squared  of  the  grid  cell  thickness 
(Eq.  ( 1 )),  resulting  in  greater  drag  in  the  bottommost  grid  cell  in  the 
80-layer  model  compared  to  the  40-layer  model. 

Average  speeds  in  the  near-bottom  layers  over  the  SE  in  the  VQS 
coordinate  model  are  small,  remaining  less  than  0.02  m  s  1  over  the 
lower  SE  and  0.01  m  s  1  over  the  upper  SE  (Fig.  7c  and  d).  In  the  VQS 
models,  both  maximum  and  mean  currents  driven  by  truncation  er¬ 
rors  stabilize  or  have  a  slight  positive  trend  after  20  days  of  integra¬ 
tion.  In  the  NVS  models,  maximum  speeds  maintain  small  positive 
trends  through  the  integration  and  do  not  tend  to  asymptote,  except 
for  the  mean  spurious  speeds  over  the  upper  SE  ( Fig.  7d)  where  mean 
currents  rapidly  reach  a  peak  value  (0.035-0.04  ms1)  and  then  le¬ 
vel  out.  These  tests  demonstrate  that  implementation  of  the  VQS 
vertical  coordinate  system  in  a  numerical  model  of  the  SE  region 
substantially  reduces  spurious  currents  arising  from  truncation  er¬ 
rors  so  that  this  type  of  model  can  reasonably  be  used  to  simulate 
the  intense  near-bottom  currents  in  the  region. 

5.  Realistic  simulation  of  deep  currents  over  the  Sigsbee 
Escarpment 

The  nested  modeling  system  is  applied  as  described  in  Section 
3.3  to  simulate  energetic  near-bottom  currents  over  the  SE  with 


characteristics  consistent  with  observations.  The  nested  modeling 
system  consists  of  the  1/120°  SE  domain  obtaining  initial  and 
boundary  conditions  from  the  outer  1/20°  CoM  domain.  First,  the 
outer  GoM  domain  is  integrated  for  10  years  and  the  output  ana¬ 
lyzed  to  identify  time  periods  of  strong,  deep  currents  over  the  SE. 

Similar  to  observations  (Donohue  et  al.,  2006),  the  GoM  model 
simulates  two  regions  with  strong  currents  over  the  SE.  At  the  east¬ 
ern  end  of  the  lower  SE,  the  maximum  speed  over  the  10-year  per¬ 
iod  of  integration  is  0.59  m  s  1  (Fig.  8a).  Over  the  upper  SE,  a  local 
maximum  of  0.38  ms  1  is  simulated  during  the  model  experiment. 
This  intensification  of  the  near-bottom  currents  over  the  SE, 
although  relatively  strong  for  such  depths,  is  somewhat  weaker 
than  has  been  observed.  One  possible  reason  is  that  the  1/20°  grid 
spacing  of  the  GoM  domain  cannot  adequately  resolve  the  small 
spatial  scales  of  the  intensified  currents  or  of  the  escarpment. 
The  nested  model  approach  refines  the  grid  to  allow  the  model 
to  more  accurately  resolve  these  large  topographic  and  velocity 
gradients  necessary  for  properly  simulating  the  strong  current 
events. 

Several  times  during  the  1 0-year  simulation,  the  SE  region  is  di¬ 
rectly  impacted  by  energetic  mesoscale  eddies.  Based  on  previous 
studies,  it  is  expected  that  stronger  near-bottom  currents  over 
the  SE  are  more  likely  to  occur  during  these  times.  One  year  of 
the  outer-domain  simulation  is  selected  for  a  one-year  nested 
model  simulation.  For  this  case  study,  relatively  strong  currents 
are  simulated  in  the  GoM  model  domain  over  the  SE  when  the  Loop 
Current  has  penetrated  far  to  the  northwest  and  separates  a  Loop 
Current  Eddy,  centered  just  southeast  of  the  SE.  During  this  year, 
several  energetic,  near-bottom,  current  events  are  evident  in  the 
GoM  model  domain  over  the  upper  SE,  with  a  maximum  current 
speed  of  0.34  ms  1  over  the  upper  SE  and  ^0.3  m  s  1  over  the  low¬ 
er  SE.  The  high-resolution  SE  model  is  initialized  from  fields  inter¬ 
polated  from  the  GoM  domain  at  the  first  day  of  the  case  study  and 
integrated  for  365  days  with  boundary  conditions  being  obtained 
from  the  GoM  domain  every  48  hours  and  output  being  saved 
every  6  hours.  A  series  of  analysis  methods  are  applied  to  the  mod¬ 
el  output  for  comparison  with  observations. 
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Fig.  6.  Maximum  speed  in  the  lower  200  m  after  20  days  of  the  test  runs  of  the 
models  with  different  vertical  discretization  (a)  80  NVS  levels;  (b)  40  NVS  levels;  (c) 
80  VQS  levels.  Speeds  shallow  than  1000  m  are  not  shown.  Isobaths  are  drawn 
every  500  m  starting  at  3500  m.  Values  of  local  maxima  are  shown. 


5. 1  Extreme  near- bottom  currents 

Compared  to  the  coarser-resolution  GoM  domain,  the  high-reso- 
lution  SE  nested  model  simulates  much  stronger  currents  along  the 
SE  (Fig.  8b).  The  intensification  of  the  near-bottom  currents  is  more 
clearly  localized  along  the  escarpment.  Locations  of  the  near-bottom 
maxima  are  similar  to  observations  (Section  2.2).  The  most  energetic 
near-bottom  (in  the  lower  250  m)  currents  over  the  upper  SE  are 
simulated  near  the  12  mooring  location  (0.81  ms  ’)  and  at  the 
southwestern  side  of  Green  Knoll  (0.98  ms1).  The  maximum 
near-bottom  current  speed  simulated  within  the  nested  domain  is 
117ms  1  and  occurs  over  the  lower  SE.  Along  the  Escarpment,  cur¬ 
rents  are  most  intense  along  the  convex  contours  of  the  isobaths 
(appearing  similar  to  capes  in  a  coastline),  while  concave  contours 
appear  to  be  “shaded  regions”  where  speeds  are  diminished.  Thus, 
current  speeds  at  two  nearby  locations  along  the  escarpment  may 
be  drastically  different  (this  feature  has  been  noticed  from  observa¬ 


tions  and  studies  cited  in  Section  2.2).  The  existence  of  such  extreme 
currents  (and  even  more  intense  currents)  at  the  upper  SE  is  sup¬ 
ported  by  observations  (Hamilton  and  Lugo-Fernandez,  2001).  At 
the  lower  SE,  the  MMS  moorings  did  not  measure  currents  stronger 
than  0.7  m  s  1  during  a  one-year  interval.  This  should  not  cast  doubt 
on  the  model  results  because  one  year  of  observations  at  the  lower 
SE  may  not  be  long  enough  to  catch  extreme  events,  which  are  very 
episodic.  During  the  MMS  observations,  no  extreme  events  occurred 
in  the  area  (Donohue  et  al.,  2006). 

5.2.  Mean  current  speeds 

Mean  deep  currents  computed  from  the  one-year  nested  model 
simulation  have  similarities  with  the  observations  collected  in  the 
region  (Section  2.2).  A  map  of  the  modeled,  mean,  near-bottom 
speed  (Fig.  8c)  shows  two  distinct  regions:  one  with  quite  weak 
(<0.1  ms  !)  flows  located  over  shallower  depths  to  the  north  of 
the  escarpment  (except  for  the  western  part  where  more  energetic 
mean  currents  are  simulated),  and  the  other  with  more  energetic 
currents  along  the  escarpment  and  over  deeper  waters.  Along  the 
Escarpment  are  a  few  localized  regions  with  mean  current  speeds 
greater  than  0.2  ms  \  One  appears  at  the  eastern  end  (around 
91  °W)  of  the  lower  SE  with  a  mean  speed  of  ^0.3  ms  \  The  upper 
SE  is  almost  quiescent  with  speeds  below  0.1  m  s  \  At  the  bottom 
of  the  slope,  the  mean  currents  are  0  13-0.18  ms  \  An  intensifica¬ 
tion  of  the  mean  currents  up  to  0.2  ms-1  is  simulated  at  the  east¬ 
ern  end  of  the  upper  SE  where  the  isobaths  converge. 

Interestingly,  the  area  to  the  southwest  of  Green  Knoll  has  very 
weak  mean  currents  (<0.1  ms  1),  which  does  not  correspond  to 
the  map  of  maximum  speeds  (Fig.  8b)  where  very  energetic  flows 
are  simulated.  Inspection  of  an  animation  of  the  near-bottom 
speeds  reveals  that,  except  for  one  very  energetic  event,  this  area 
is  shaded  by  the  knoll,  which  obstructs  the  predominant  flow  along 
the  escarpment  from  the  east.  Near  the  12  mooring  location,  strong 
flows  are  more  frequent,  resulting  in  higher  mean  currents 
(although  the  simulated  maximum  speed  is  smaller  at  12  than  at 
Green  Knoll).  Intensification  of  the  currents  to  the  northwest  and 
southeast  of  Green  Knoll  indicates  flow  intensifying  as  it  is  directed 
around  the  topographic  feature. 

5.3.  Mean  velocity  vectors  and  standard  deviation  ellipses 

Velocity  vector  means  and  standard  deviations  are  computed  for 
select  grid  points  over  the  SE  region  (Fig.  9).  Along  the  lower  and 
upper  SE,  the  mean  current  has  a  well-pronounced  westward  com¬ 
ponent.  The  standard  deviations  (STD)  of  the  mean  vectors,  dis¬ 
played  as  ellipses  at  the  heads  of  the  vectors,  characterize  the 
variability  of  the  mean  currents  at  the  given  locations.  Variability 
of  the  simulated  mean  current  vectors  agrees  well  with  characteris¬ 
tics  of  the  measured  currents  (discussed  in  Section  2.2  and  Donohue 
et  al.  (2006)).  The  prominent  feature  of  the  STD  ellipses  in  the  model 
is  the  along-isobath  orientation  of  the  major  axes  along  the  SE  where 
the  isobaths  converge,  which  is  indicative  of  quasi-rectilinear  behav¬ 
ior  of  the  currents  over  the  steep  slope.  Over  the  flatter  regions,  the 
ellipses  are  more  circular  with  their  major  axes  rotated  at  a  larger  an¬ 
gle  with  respect  to  the  isobaths.  The  semi-major  axes  of  the  ellipses 
are  comparable  in  magnitude  to  the  mean  vectors.  This  shows  that 
the  simulated  currents  at  the  upper  and  lower  SE  are  highly  variable. 
The  velocities  range  from  near-zero  to  very  strong  currents  directed 
westward  along  the  isobaths.  The  STDs  of  the  model  vector  velocities 
suggest  that  flows  rarely  reverse  in  this  area. 

5.4.  Vector  velocity  time  series 

Time  series  of  the  simulated  currents  from  the  high-resolution 
nested  model  domain  (Fig.  10)  are  analyzed  for  selected  points 
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Fig.  7.  Time  series  of  the  maximum  (upper  row)  and  spatially-averaged  (bottom  row)  speeds  over  the  Lower  (left  column)  and  Upper  (right)  SE  for  test  model  runs  with  3 
different  vertical  grids;  80  non-vanishing  a  layers  ("80  NVS"),  40  non-vanishing  a  layers  (“40  NVS"),  and  80  vanishing  (quasi)  a  layers  ("80  VQS").  The  speeds  are  for  the  near¬ 
bottom  coordinate  surface. 


across  the  lower  and  upper  SE.  The  points  are  chosen  in  the  middle 
of  the  escarpment  along  90°W  at  the  upper  SE  and  91°W  at  the 
lower  SE.  These  points  are  located  in  the  areas  with  high  maximum 
near-bottom  speeds  (see  Fig.  8b).  For  each  grid  point,  vector  time 
series  are  plotted  showing  currents  at  90  m  below  the  surface 
and  approximately  800,  250,  and  50  m  above  the  bottom.  These 
time  series  have  many  similarities  to  velocity  time  series  collected 
during  the  MMS  observational  study  discussed  in  Section  2.2: 

•  In  the  model,  the  near-surface  and  lower-layer  vectors  show  no 
relationship,  meaning  that  the  upper-  and  lower-level  flows  are 
not  coherent.  Currents  are  coherent  in  the  lower  water  column. 

•  Currents  are  intensified  towards  the  bottom. 

•  Due  to  the  frictional  bottom  boundary  layer,  the  maximum  flow 
in  the  lower  water  column  is  separated  some  distance  above  the 
bottom.  In  the  model  experiment,  the  maximum  of  the  current 
speed  during  the  most  energetic  events  occurs  at  approximately 
200-400  m  above  the  bottom. 

A  noticeable  difference  between  the  simulated  deep  currents 
and  the  observations  at  the  12  mooring  location  is  a  more  persis¬ 
tent  westward  orientation  of  the  velocity  vector  in  the  model.  Re¬ 
corded  current  vectors  exhibit  rotational  behavior  and  more 
frequently  change  direction  from  westward  to  eastward  (compare 
Fig.  10b  with  Fig.  2  in  Hamilton  and  Lugo-Fernandez,  2001).  Unfor¬ 
tunately,  at  this  stage  it  is  impossible  to  state  whether  this  differ¬ 
ence  is  related  to  artifacts  of  the  simulation  or  differences  in  the 
dynamical  regimes  from  this  one-year  simulation  and  the  time 
period  of  the  observations. 

Time  vs.  depth  diagrams  of  speed  (Fig.  1 1 )  more  clearly  demon¬ 
strate  features  noted  from  the  velocity  time  series.  The  diagrams 
show  that  extreme  events  are  episodic  and  not  regular  at  the  SE. 
The  lower  SE  (Fig.  1  la)  has  more  strong  events  over  the  simulation 
(during  days  80  through  180  of  the  integration  and  a  single  event 
around  day  310).  The  upper  SE  (panel  b),  has  a  single  strong  event 
(occurring  over  a  short  time  period  centered  near  day  125)  and 
several  less  energetic  episodes.  Both  at  the  lower  and  upper  SE, 
the  deep  currents  are  decoupled  from  the  upper  ocean.  Intensifica¬ 
tion  of  the  near-bottom  currents  is  typically  seen  with  only  modest 
speeds  at  the  surface.  The  time  vs.  depth  diagrams  of  speeds  of  the 


simulated  currents  vividly  demonstrate  near-bottom  intensifica¬ 
tion  of  the  flow  and  the  existence  of  shear  due  to  the  bottom 
boundary  layer.  Along  both  the  lower  and  upper  SE,  the  current 
is  intensified  a  few  hundred  meters  above  the  bottom  and  then 
slowly  decays  towards  the  upper  ocean  (above  the  main  ocean 
thermocline,  very  intense  currents  are  associated  with  upper- 
ocean  circulation  features). 

The  simulated  currents  at  the  upper  SE  have  fairly  good  resem¬ 
blance  to  observed  oscillations  at  the  12  mooring  during  the  strong 
event  (compare  panels  a  and  b  in  Fig.  12).  Both  time  series  show 
shear  near  the  bottom,  likely  due  to  the  bottom  boundary  layer, 
with  flow  intensification  at  200-400  m  above  the  bottom.  The 
dominant  period  of  the  observed  oscillations  is  approximately 
10  days,  which  matches  noticeably  well  with  the  model  results. 
In  the  simulated  time  series,  the  10-day  oscillation  seems  to  be 
superimposed  over  a  slower  mode. 

5.5.  Wavelet  analysis 

Analyses  of  the  observations  reveal  that  the  dominant  periodic¬ 
ity  of  the  low-frequency  signal  in  the  deep-ocean  currents  over  the 
SE  is  10-40  days  (Donohue  et  al.,  2006).  The  time  series  of  the  sim¬ 
ulated  currents  are  analyzed  to  estimate  the  time  scale  of  the  var¬ 
iability  of  the  deep-ocean  dynamics  in  the  model.  The  velocity  time 
series  from  the  one-year  simulation  are  nonstationary.  The  oscilla¬ 
tions  do  not  reveal  a  robust  periodicity  throughout  the  record. 
Thus,  applying  a  wavelet  transform  to  the  time  series  aids  in  ana¬ 
lyzing  periodic  behavior  of  the  variability  of  the  near-bottom  cur¬ 
rents  over  the  escarpment. 

The  along-isobath  velocity  component  is  transformed  using  the 
Morlet  wavelet  (Goupillaud  et  al.,  1984): 

g(t)  =  cmn  l/4e  'J/2(e'"’'  -  e  ,"‘/2)  (5) 

where  ca)  is  the  normalization  constant: 

C,„  =  (l  +e  '"2  -2e  w'4)  1/2.  (6) 

To  reduce  ringing,  each  end  of  the  time  series  is  extended  by 
adding  a  trigonometric  taper  (1  sin<£)  where  <f>  is  chosen  such 
that  the  taper  vanishes  over  the  extended  interval  (100  days). 
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Fig.  8.  (a)  Maximum  near-bottom  (lower  250  m)  speed  over  the  SE  region  during 
10  years  of  the  GoM  model  run.  The  arrows  indicate  locations  of  local  speed 
maxima  over  the  10-year  integration  at  the  lower  and  upper  SE.  The  grey  area  is 
where  the  seafloor  is  shallower  than  1000  m.  (b)  Maximum  near-bottom  (lower 
250  m)  speed  over  1  year  of  the  SE  model  run.  (c)  Near-bottom  speed  averaged  over 
1  year  of  the  integration. 

Over  the  lower  SE,  the  wavelet  transform  indicates  oscillations  of 
a  more  complex  nature  (Fig.  13a)  than  those  at  the  upper  SE 
(Fig.  13b).  Over  the  deeper  SE,  the  oscillations  are  less  regular  and 
persistent  with  frequency  rapidly  shifting  higher  and  lower.  Prior 
to  model  year-day  120,  the  site  is  influenced  by  two  wave  patterns: 
one  with  low-frequency  motions  with  periods  of  70-90  days,  and 
the  other  with  higher-frequency  oscillations  of  approximately  20- 
day  periodicity.  These  faster  oscillations  are  more  prominent  at 


the  lower  SE  than  at  the  upper  SE.  After  about  1 30  days  of  the  nested 
model  run,  the  frequencies  of  the  fast  and  slow  modes  begin  to  con¬ 
verge  and  the  dominant  periodicity  gradually  shifts  to  35-65  days 
by  day  200.  During  days  210-270,  there  are  two  prominent, 
although  weakened,  modes  of  oscillations  with  more-defined  peri¬ 
odicity  (the  maxima  are  more  localized).  The  slower  mode  shows 
45-60  day  periodicity,  dissipating  toward  the  end  of  the  integration, 
and  faster  oscillations  of  about  20-day  periodicity  are  well  identified 
in  the  diagram. 

For  the  upper  SE,  the  wavelet  diagram  (Fig.  13b)  reveals  a  dom¬ 
inant  wave  pattern  with  the  period  of  oscillation  ranging  from  35 
to  65  days.  This  pattern  is  surprisingly  distinct  and  robust  through¬ 
out  the  simulation.  The  frequency  of  these  oscillations  only  slightly 
drifts  lower  during  the  simulation.  The  maxima  of  the  wavelet 
transform  coefficients  are  well  localized  in  this  frequency  interval, 
meaning  that  these  oscillations  are  regular  and  should  be  easily 
identifiable  by  spectral  analysis.  There  are  also  a  clear,  slowly- 
oscillating  mode  (100-120  days)  and  weak  oscillations  with  IQ- 
30  day  period.  The  slowly-oscillating  pattern  may  be  associated 
with  slow  changes  in  the  character  of  the  upper-ocean  mesoscale 
dynamics.  The  period  of  the  faster  modes  (~10  and  30  days)  agrees 
with  the  observed  periodicity  in  the  area  (Donohue  et  al.,  2006),  It 
should  be  noted  that  this  faster  mode  intensifies  during  the  time 
period  that  intense  currents  are  simulated  in  the  area  near  model 
year  days  110-130.  Similar  10-day  oscillations  are  observed  in 
the  time  series  of  the  speed  of  the  recorded  currents  at  the  12 
mooring  (Fig.  12b). 

The  wavelet  analysis  of  the  nonstationary  velocity  time  series 
identifies  differences  in  the  character  of  the  oscillations  between 
the  upper  and  lower  SE,  and  identifies  multiple  modes  of  oscilla¬ 
tion  The  periodicities  of  the  oscillations  tend  to  shift  during  the 
year-long  integration,  but  the  relatively  fast  modes  are  evident 
throughout.  The  enhancement  of  the  amplitude  of  the  fast  oscilla¬ 
tions  during  the  most  energetic  events  is  consistent  with  previous 
studies  and  adds  to  the  confidence  that  the  model  is  indeed  repro¬ 
ducing  the  observed  dynamics  over  the  escarpment. 


6.  Discussion  and  conclusions 

This  modeling  effort  was  motivated  by  a  series  of  observations 
of  extreme  deep  currents  over  the  SE.  The  physical  phenomenon 
responsible  for  the  intense  near-bottom  flows  has  not  yet  been 
fully  explained  and  requires  a  comprehensive  (observational,  ana¬ 
lytical,  and  numerical)  investigation  for  the  benefit  of  both  scien¬ 
tific  and  engineering  applications.  This  study  is  the  first  known 
attempt  to  simulate  extreme  near-bottom  currents  over  the  SE 
with  a  numerical  model.  A  new  VQS  vertical  grid  has  been  devel¬ 
oped  to  overcome  difficulties  faced  by  traditional  ocean  models 
for  studying  deep  currents  over  a  region  with  such  complicated 
topography  and  steep  slopes.  Model  test  experiments  have  demon¬ 
strated  that  the  new  vertical  grid  results  in  smaller  numeric  errors 
for  the  simulation  over  the  steeply  sloping  bottom  compared  to  a 
traditional  sigma-layer  model.  Implementation  of  the  new  VQS 
vertical  grid  with  up  to  80  unevenly  stretched  vertical  layers  pro¬ 
vides  accurate  representation  of  the  SE  slope  and  the  highest 
(known  to  date)  vertical  resolution  over  the  escarpment. 

A  nested  modeling  system  has  been  developed  for  this  study  to 
adequately  resolve  the  small  spatial  scales  of  the  escarpment  and 
the  observed  intense  currents.  With  this  approach,  a  very-high-res- 
olution  (approximately  800-900  m )  model  domain  can  be  applied  to 
the  SE  region,  while  still  being  influenced  by  realistic,  upper-ocean, 
mesoscale  circulation  features  of  the  Gulf  of  Mexico  simulated  by 
the  outer  nest.  The  dynamic  fields  simulated  using  this  nested  model 
approach  have  properties  that  agree  well  with  the  prominent  fea¬ 
tures  of  the  observed  currents  over  the  Escarpment  (Sections  2.3): 
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Fig.  9.  Near-bottom  mean  velocity  vectors  and  standard  deviation  ellipses  over  the  SE  estimated  from  1  year  of  data  from  the  Sigsbee  Escarpment  model  Locations  of  the 
mean  vector  estimates  are  marked  with  asterisks.  Grey  isobaths  are  drawn  every  200  m.  and  the  1500  and  2000  m  isobaths  are  shown  by  black  contours. 
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Fig.  10.  Time  series  of  daily  velocity  vectors  for  one  year  of  the  SE  model  run.  The  column  (a)  and  (b)  show  velocity  at  grid  points  at  the  lower  and  upper  SE.  respectively  (the 
coordinates  are  in  the  heading).  The  bottom  depth  and  coordinates  of  the  grid  points  are  given  in  the  headings.  The  horizontal  axis  is  the  model  time  in  days.  For  each  point, 
time  series  are  presented  at  4  depth  levels  indicated  in  the  diagrams  as  "h".  Note  the  orientation  of  the  north  direction  shown  in  the  upper  panels. 


Mean  currents.  Along  the  SE,  the  simulated,  near-bottom,  mean 
currents  are  directed  predominantly  southwestward  with  major 
axes  of  the  STD  ellipses  oriented  at  small  angles  to  the  isobaths, 
indicating  quasi-rectilinear  oscillations.  Speeds  of  the  mean  cur¬ 
rents  in  the  model  range  from  0.1  to  0.3  m  s  \  Currents  are  weaker 
over  the  shallower  topography  to  the  north  of  the  SE  (less  than 
0.1  ms  ’).  Along  the  slope,  intensification  of  the  mean  currents 


is  observed  at  locations  seemingly  influenced  by  topographic  irreg¬ 
ularities  along  the  lower  and  upper  SE. 

Extreme  currents.  Several  events  with  intense  near-bottom  cur¬ 
rents  are  simulated  along  the  SE.  The  strongest  currents  are  located 
near  locations  where  similar  extreme  currents  have  been  observed. 
At  the  upper  SE,  the  strongest  simulated  speed  is  0.81  ms  ’.Stron¬ 
ger  speeds  up  to  1.17  m  s  1  are  simulated  along  the  lower  SE. 
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Fig.  11.  Time  (horizontal  axis)  vs  depth  (vertical  axis)  diagrams  of  speeds  of  the  simulated  velocities  at  the  lower  (a)  and  upper  (b)  SE.  Coordinates  of  the  locations  are 
indicated  in  the  headings. 
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Fig.  12.  Current  speed  at  different  depth  levels  in  the  model  (a)  and  12  mooring  observations  (b)  at  the  upper  SE.  The  model  time  series  in  (a)  is  48  HLP  filtered  and  the  grid 
point  is  located  close  to  the  12  mooring  (90°W,  27.2 1°N.  same  as  Fig.  1 1b).  (b)  Current  speed  of  36  HLP  filtered  velocity  data  from  the  12  mooring.  The  arrows  indicate  peaks  of 
a  quasi- 10-day  period  signal. 


Periodicities  of  oscillations.  From  inspection  of  the  time  series 
of  current  speed  and  the  wavelet  analysis  of  the  model  velocity, 
the  estimated  periodicity  of  the  low-frequency  signal  during  strong 
events  in  the  region  is  10-30  days  in  agreement  with  observations. 
The  slower  modes  of  the  simulated  currents  with  periods  40- 
60  days  observed  in  the  wavelet  diagrams  can  also  be  seen  in  some 
of  the  spectra  from  the  Donohue  et  al.  (2006)  study.  Analysis  of  the 
simulation  suggests  that  these  slow-oscillating  modes  are  not 
attributed  to  the  strong  events. 

The  model  development  and  evaluation  described  in  this  work 
is  intended  to  produce  a  modeling  methodology  that  can  be  ap¬ 
plied  to  future  studies  of  the  dynamical  aspects  these  extreme 
events  over  the  SE.  Nevertheless  some  preliminary  results  related 
to  the  dynamics  of  the  deep  currents  over  the  SE  from  analyses 
of  the  model  experiments  can  be  mentioned  here. 

•  The  modeling  results  show  dominantly  quiescent  conditions 
over  the  area,  but  with  occasional  strong  events  when  an  eddy 
migrates  over  the  region,  especially  over  the  upper  SE.  This 
result  agrees  with  observations  that  suggest  a  linkage  between 
the  near-bottom  extreme  events  and  intense  eddy  activity  in 
the  vicinity.  However,  some  strong  events  are  simulated  over 


the  lower  SE  away  from  the  upper  ocean  eddy  activity.  Thus, 
the  model  results  suggest  the  possibility  for  different  genesis 
of  the  extreme  events  over  the  lower  (local  and  remote  genesis) 
and  upper  (local)  SE. 

•  Strong  events  simulated  over  the  SE  appear  to  propagate.  It  is  not 
clear  yet  whether  the  propagation  of  the  strong  currents  along  the 
slope  is  a  manifestation  of  a  topographically-trapped  wave  or  a 
deep  eddy  impinging  the  slope.  Both  possibilities  are  supported 
by  observations  (Donohue  et  al.,  2006).  The  applicability  of 
Rhine’s  (1970)  TRWs  theory  to  the  studied  case  should  be  revis¬ 
ited.  The  slope  of  the  SE  is  steeper  than  slopes  permitted  in  Rhines 
(1970)  theory  for  TRW-like  motion.  Also,  bottom  friction  is  not 
taken  into  account  in  this  theoretical  description  of  TRWs. 

•  Directly  over  the  SE,  flows  in  the  deep  ocean  are  predominantly 
not  coherent  with  flows  at  the  surface. 

•  Simulated  strong  events  are  localized  (in  agreement  with  obser¬ 
vations,  e.g.,  Hamilton,  2007),  suggesting  that  there  are  certain 
areas  along  the  escarpment  favored  for  flow  intensification. 

•  The  maximum  speed  of  the  modeled  flows  is  achieved  at  some 
height  (~100-200m)  above  the  bottom,  which  suggests  the 
influence  of  the  bottom  boundary  layer  on  the  vertical  structure 
of  the  extreme  currents. 
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Fig.  13.  Morlet  wavelet  transform  coefficients  (amplitude)  of  the  time  series  of  the 
along-isobath  component  of  the  near-bottom  (2S0  m  above  the  bottom)  velocity  at 
the  lower  (a)  and  upper  (b)  SE.  The  location  of  the  grid  point  at  the  upper  SE  is 
27.2N,  90W,  and  at  the  lower  SE  is  26. IN.  91 W.  The  mesh  indicates  the  cone  of 
influence  at  the  endpoints,  meaning  the  wavelet  coefficients  are  not  significant  in 
that  area.  The  colors  represent  values  of  the  coefficients  relative  to  the  minimum 
and  maximum.  The  vertical  axis  is  the  scale  length  of  the  wavelet  converted  into 
periods.  The  horizontal  axis  is  the  translation  parameter  describing  the  position  of 
the  wavelet  on  the  time  axis. 


Although  the  results  of  the  simulation  demonstrate  that  the 
newly-developed  model  has  met  the  goals  of  the  study,  the  analy¬ 
sis  of  the  numerical  experiments  highlights  aspects  of  the  model 
that  potentially  require  further  investigation  and  might  yet  be  im¬ 
proved  to  achieve  better  simulations  of  the  deep-ocean  dynamics 
and  bottom-intensified  features  over  the  SE.  Of  concern  is  the  role 
that  the  ocean  bottom  boundary  layer  seems  to  play  in  altering  the 
vertical  structure  of  the  bottom-intensified  currents.  More  atten¬ 
tion  should  be  placed  on  accurately  modeling  the  boundary  layer 
by  investigating  the  consequences  of  the  turbulence  closure 
scheme  and  the  parameterization  of  bottom  roughness.  Changing 
parameters  of  the  bottom  friction  drag  coefficient  (1)  can  affect 
the  location  of  the  maximum  relative  to  the  local  bottom,  though 
it  is  unlikely  that  it  will  drastically  change  the  intensity  of  the 
current. 

It  is  possible  that  the  horizontal  resolution  of  the  SE  regional 
model  is  still  inadequate  for  accurately  representing  the  steep 
topographic  features  and  this  should  be  addressed.  The  effect  of 
the  vertical  resolution  in  the  model  on  simulating  the  strong 
near-bottom  current  has  not  been  thoroughly  investigated.  Finally, 
a  hindcast  simulation  of  a  strong  event  needs  to  be  performed  to 
more  directly  compare  the  model  results  with  the  observations 
and  validate  the  model.  Nevertheless,  this  modeling  methodology 
can  be  applied  to  significantly  advance  understanding  the  dynam¬ 
ics  over  the  Sigsbee  Escarpment  region. 
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